Stability of charge inversion, Thomson problem and application to electrophoresis 
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We analyse charge inversion in colloidal systems at zero temperature using stability concepts, and 
connect this to the classical Thomson problem of arranging electrons on sphere. We show that for a 
finite microion charge, the globally stable, lowest energy state of the complex formed by the colloid 
and the oppositely charged microions is always overcharged. This effect disappears in the continuous 
limit. Additionally, a layer of at least twice as many microions as required for charge neutrality 
is always locally stable. In an applied external electric field the stability of the microion cloud is 
reduced. Finally, this approach is applied to a system of two colloids at low but finite temperature. 

PACS numbers: 82.70.Dd, 41.20.-q 
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I. INTRODUCTION 

We study analytically and numerically charged col- 
loidal particles in the presence of electrolyte solutions in 
the low-temperature limit. This is equivalent to strong 
electrostatic coupling at finite temperature jjj. Already 
at room temperature a colloidal system thus is in the low 
temperature regime, provided the electrolyte consists of 
multivalent ions. In this limit, a certain number of ions 
condenses onto the surface of the colloid. For consistency, 
we will call the charged colloids macroions, and the ions 
of the electrolyte microions. For the general macroion 
problems, this behaviour has recently attracted a lot of 
attention especially due to its importance in biological 
systems § | g 1 |J, ||. 

Macroion complexes exhibit numerous counter- 
intuitive phenomena. The most pronounced one is 
referred to as overcharging or charge inversion. A cer- 
tain number of microions is needed to condense for the 
entire complex (macroion and condensed microions) to 
become charge neutral. Sometimes an excessive number 
of microions condenses [g g, |, g ||, ||[ || ||, [jj, 
g 0, and the complex, the "dressed" 

macroion, acquires an effective charge that is opposite 
in sign to that of the bare macroion. This phenomenon 
is beyond the standard Debye-Hiickel [|2|, |23|, Q 
and Derjaguin-Landau- Verwey-Overbeek (DLVO) 
theories p5[ . 

Another phenomenon is auto-ionisation jl6|, |l7|, |l9| . 
This means that one macroion transfers some of its mi- 
croions to another macroion, so that the first one becomes 
undercharged while at the same time the second one be- 
comes overcharged. Also analysed in the literature is the 
important question of transport in an external electric 
field. A bare macroion will move in the direction deter- 
mined by its own charge. The binding of microions to it 
can under certain conditions reverse the direction. 

These results have mostly been arrived at by molecu- 
lar dynamics (MD) and Monte Carlo (MC) simulations. 
Simulations offer the advantage that finite temperature 
can be taken into account in a natural way. Purely an- 
alytical approaches have to resort to relatively compli- 



cated starting points since mean-field theories are insuf- 
ficient ||^, gij. A successful and often used approach 
is the model of a two-dimensional Wigner crystal which 
becomes exact at zero temperature and very large num- 
ber of microions per macroion ]3^ . 

At zero temperature — without any restrictions on the 
number of microions — the problem is, however, directly 
related to the classical "Thomson problem" of finding 
stable configurations of N mutually repelling electrons 



!7j]. While the original 
'udding" model for the 



on the surface of a sphere J26| , 
problem was about the "Plum 
atom, where N particles are confined inside a homoge- 
neously charged sphere, both problems actually arc iden- 
tical since the repelling interaction will push all particles 
inside the sphere toward its surface. This similarity be- 
tween overcharging and the Thomson problem seems to 
have been largely unnoticed in the macroion literature. 

In this paper we will make use of the Thomson prob- 
lem to derive rigorous bounds for the phenomena dis- 
cussed above. At zero temperature, our results are exact 
and hence an improvement on previously known results 
from Wigner crystal theory. Compared to the results of 
simulations, our derivation suffers from our inability to 
include finite temperature in an exact way. However, our 
method outperforms previous ones in both the ease of 
the method (both conceptually and numerically) and in 
allowing to treat many phenomena in a single consistent 
way. 

This paper is organised as follows. In Sec. |H| we 
summarise the so-called primitive model. This model 
is used in basically all studies of macroion complexes. 
In Sees. |Tn|-|v| we analyse the stability of overcharged 
macroion complexes. We use concepts from dynamical- 
systems theory to show that two different stability prop- 
erties exist, global stability (Sec. Ill) and local stability 
(Sec. IV). In Sec. VI we move on to the question of a 
macroion in an applied external field, i.e., electrophore- 
sis. We will discuss the auto-ionisation of macroions in 



Sec. VII. We conclude in Sec. VIII 
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II. MODEL 

We consider a spherical macroion of charge Q and ra- 
dius R mac , surrounded by N spherical microions of charge 
q and radius i? m i C . Q and q are of opposite signs, and in 
the following we assume Q < 0. The macroion is fixed 
at the origin and the N microions are distributed at po- 
sitions r\, i — 1, . . . , N. The total electrostatic energy V 
for a particular configuration is then given by 
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where the first sum accounts for the mutual repulsion of 
the microions, and the second for the attraction between 
the macroion and each microion. The effect of solvent is 
included through an effective dielectric constant e. Short- 
range pairwise repulsion is taken into account by hard- 
core interaction 
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where v(r) — + oo for r < and zero otherwise. 
Equations (|l|) and (|J) comprise the so-called primitive 
model 01. 



III. GLOBAL STABILITY 

We approach the pro blem by using the well-known 
Earnshaw's theorem [Eg which states that there can be 
no stable state in a system with only electrostatic inter- 
actions present. For stable configurations to exist, short- 
range repulsive forces must be present in addition to the 
long-range Coulomb ones. For our system, the short- 
range forces are due to hard-core interaction [Eq. (0)]. 
Earnshaw's theorem thus restricts stable configurations 
to have all microions at a distance R = i? m ac + i? m ic away 
from the centre of the macroion. 

The condition = R allows us to simplify Eq. (Q) to 
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with the normalised coordinates 



r' ' 



Equation (||) 

no longer describes the energy of an arbitrary arrange- 
ment {fi} of particles but the energy of any stable ar- 
rangement instead. Next, we introduce the function 



f(N) 
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with 



(4) 



where we demand that the coordinates are the ones 
for the lowest energy state with N microions around the 
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FIG. 1: Comparison of Eq. (^) [solid line] with the result of 
a numerical computation of Eq. (ft) [open circles] . 



macroion. Thus, the coordinates are completely defined 
by N. 

At the ground state Eq. (Q) becomes minimised. The 
complete solution can be computed numerically very ef- 
ficiently ^9|. Furthermore, the functional form of f(N) 
is known to excellent precision |27j to be 

N 2 

f(N) = — - cN 3/2 with c = 0.5510 . (5) 

This formula is easy to understand when one notices that 
the first term is the energy of a continuous layer of charge 
on a sphere of unit radius while the second term is the 
self-energy correction due to discrete microions which can 
be shown to be proportional to TV 3 / 2 . A comparison of 
Eq. (H) and a numerical solution of the exact formula 
[Eq. (H)] is shown in Fig. |l|. 

By using the condition |f | = R for all i we have com- 
pletely accounted for the hard-core interaction between 
the macroion and the microions. We can neglect the 
hard-core interaction between microions since they repel 
each other as their charges have the same sign. Collecting 
results, the potential energy of the lowest energy state for 
a macroion surrounded by N microions is given by 



V(N) = 
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Earnshaw's theorem gives a necessary but not suffi- 
cient criterion for the stability of a system. Furthermore, 
it states that an unstable microion is immediately pushed 
to infinity. Let us consider a macroion and N microions 
where we place the microions at arbitrary positions - 
not necessarily on the macroion. Due to Earnshaw's 
theorem, M of them will attach to the macroion while 
k = N — M will escape to infinity. All N will go to the 
macroion if it is the state lowest in energy, i. e., if 



V(N) < V(M) V < M < N 



(7) 



This condition is much stronger than the simple condition 
V(N) < 0, since the latter only prevents all microions 
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FIG. 2: Globally stable overcharging Q g i b as function of the 
ratio of the charges of macro- and microion, computed from 
Eq. (S). 



from escaping simultaneously while Eq. (Q) also prevents 
some from escaping. 

Due to Earnshaw's theorem, all stable solutions are 
enumerated by the number of microions, and we simply 
have to find the number yielding the lowest energy. Since 
Eq. @ possesses only a single extremum for given pa- 
rameters q, Q and R, we can simply use dV(N)/dN = 0. 
That yields 



iV glo b = ^ + c c + — 




7^i— > (8) 



where we have used the assumption Q < introduced 
above. 

The first term \Q\/q gives the naive result that a com- 
plex consisting of macroion and layer of microions should 
be charge neutral. The other two terms give the excess 
bound microions. The maximum stable overcharging in 
terms of charge is Q g i ob = qN gloh - \Q\, see Fig. |. 

The stability criterion that we derived applies to an 
arbitrary initial placement of microions. For this reason, 
this kind of stability is referred to as global stability, see 
any textbook on nonlinear dynamics, e. g. Ref. pO. 



IV. LOCAL STABILITY 

In addition to global stability, there exists the con- 
cept of local stability. While global stability states that 
the microions will move to the macroion independent of 
their initial positions, local stability means that they will 
stay at the macroion if they have initially been placed 
there. The system is locally stable (but not globally) if 
the system could lower its energy by transferring one (or 
more) microions from the macroion to infinity but in do- 
ing so would need to cross an energy barrier. Since we 
are restricting ourselves to classical physics at zero tem- 
perature, it is impossible to cross such a barrier and the 



microions would stay touching the macroion forever — if 
prepared with this initial condition. 

To calculate the condition for the existence of such 
a barrier we move particle k slightly away from the 
macroion by a distance A, keeping all other microions 
on the surface of the macroion. If this move increases 
the potential energy, resulting in a restoring force, the 
system is locally stable. 

We label all quantities in the perturbed state by a 
prime, hence r' k = R + A and r[ = r; = R V I ^ k, 
and we introduce the abbreviation dki = |r*fc — n|- Since 
we need to consider only small A, we can use a series 
expansion, with the result 

1 1 1 * ,1 1 1 * 



d' kl dki 2Rdki 



Inserting this into Eq. (JTJ) gives 



V = V- 



q 



AireR 



I...N 

^ 2d kl T R 



A = V--^A . (10) 



The system is locally stable if and only if the expression 
in brackets is negative for all k since then an increase in 
A will increase the potential energy. Thus, the condi- 
tion for local stability is h k < V k. Since the lowest 
energy arrangement of the particles rarely is completely 
symmetric, this yields the necessary but not sufficient 
condition {hk)k < 0, where this average is over all pos- 
sible particles k. Noting that J2k hk = V(N), this gives 
the necessary condition for local stability 



V(N) < 



fll) 



It should be noted that this simple form for the condition 
is a coincidence, and for other systems V(N) < has not 
necessarily relation to local stability. With the help of 
Eq. (^) this condition can be converted into an upper 
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FIG. 3: Number of locally stable microions as a function of 
the ratio of the charges. The solid line is the analytic upper 
bound Eq. (^2|), while the dashed line is a numerical solution 
of the exact formula Eq. (|lo|). 
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bound for the number of microions that can be bound 
locally stable, 



AT loc = 2^ + 2c 




(12) 



In terms of charge this is Qi oc = qN\ oc — \Q\. 

To check the difference between Eq. ( p"2| ) and the ex- 
act solution, we have numerically computed the lowest 
energy state as a function of N, and from that deter- 
mined the largest hk for each N. The result in Fig. || 
shows that hardly any difference between the two values 
can be seen. This does not come as a surprise since the 
differences between hk for different k are small as the re- 
pelling forces among the microions try to make all mutual 
distances as equal as possible. 



V. SUMMARY OF STABILITY CONCEPTS 

For finite q, the number of microions that are bound 
globally stable is always larger than the value N — \Q\/q, 
i. e., the macroion is overcharged. In the continuous limit 
q — > this effect disappears. In contrast, the number 
of locally stable bound microions is at least twice the 
amount needed for charge neutrality, and this effect per- 
sists even in the continuous limit. Fig. || shows the dif- 
ferent regimes as a function of the charges of macroion 
and microions. 

We now want to put our results into perspective of 
previous results on overcharging [|[ |[ ||, ^, [l(], [ll], [l^, 



U 14 



H If, U, H, §(J |2J. At finite T, there are 
only few simulations done in the geometry employed in 
this paper, e. g., Ref. [2l]. Most analytical work focuses on 
T = 0, both due to simplicity and allowing one to focus 
on the influence of correlations. We employ the same 
approach. 
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FIG. 4: Number of microions that can be bound globally 
stable or locally stable to a macroion. Nq/\Q\ = 1 is the 
number of microions expected from charge neutrality. 



The main advantage of our approach is the adoption 
of the exact Thomson model as source for the potential 
energy V(N) whereas previous papers used an energy 
estimate for V(N) derived from strongly-correlated liquid 
and Wigncr crystal theory |, |TJ. The Wi gner crystal 
theory contains a parameter a whose value can either 
be determined from simulations as a function of N - 
which trivially leads to a self-consistent result — or needs 
to be fixed with the analytical value for a known for a 
two-dimensional Wigner crystal. In the latter case, this 
results in an error of up to order 10% for the computed 
energies |l6| . In contrast the error of Eq. (|J) is negligible 
(less than 10" 4 ). 

We apply methods from nonlinear dynamics, using the 
concepts of local and global stability. The existence of 
these two different stability properties seems to be unno- 
ticed in macroion literature. For example, the criterion 
by Messina et al. agrees with our global stability cri- 
terion up to the differences caused by their choice for 
V(N). The concept of local stability, however, is also 
an important one as can be seen, for example, in the 
electrophoresis setup treated in the following section. 

Finally, for numerical calculations we employ a min- 
imisation scheme. Only a few different initial conditions 
are necessary to make sure that the algorithm does not 
become stuck in a local minimum. This is in contrast to 
MD simulations which suffer from the slowing down of 
the dynamics at low temperature. 



VI. MACROION IN AN EXTERNAL ELECTRIC 
FIELD 

Let us consider a macroion with N microions in an 
external field £. We will restrict ourselves to the case of 
homogeneous external field so that the dipole and higher 
moments of the macroion complex are irrelevant. Typical 
electrophoresis experiments are done in the presence of a 
homogeneous field. 

The total force acting on the complex becomes simply 



F=(qN + Q)S 



(13) 



The complex will thus move in the same direction as the 
bare macroion if the macroion is undercharged, it will 
move in the opposite direction if it is overcharged, and it 
will remain at rest if it is charge neutral. 

Having a macroion complex with given N, q and Q, 
the interesting question is not in which direction the ef- 
fective force acts [since that question is trivially answered 
by Eq. ([jjj)] but rather whether the forces become so 
large that the system disintegrates. (For a system with 
only gravity, this problem is referred to as stability under 
tidal forces.) This was noted earlier when the dependence 
of the mobility of the macroion complex on an applied 
field was analysed |TJ, |3l| . While it was shown for a 
few examples that some microions are "ripped off" the 
macroion, no systematic study of the stability criterion 
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under an applied external field has been done (to the 
author's knowledge). Here, we aim to fill this gap. 

With an applied external field £, Eq. (|j) has to be 
extended to 
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with the reduced electric field 
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(15) 



The concept of globally stability introduced in Sec. 1 
cannot be applied in the presence of an external field 
since the potential energy is not bounded from below [j33| , 
and local stability of the complex is the relevant concept. 
Again, we move particle k by a distance A <C R away 
from the macroion. To first order the potential energy 
V' of the new state then becomes 



V' = V- 



AireR 2 
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(16b) 



The third term in the bracket of Eq. ( |l6a|) is the difference 
to Eq. ( [To| ) . It describes the interaction with the external 
field and depends on the angle between the position of 
the particle and the external field. 

As in Sec. IV the macroion complex is locally stable 
if an d on ly if hk < for all k. A closer inspection of 
Eq. (16a) and comparison to Eq. (14) shows that the big 
bracket no longer is directly related to the energy of the 
fc-th particle (as it was in Sec. IV) as the sign in front of 
S-f% is inverted @. 

We have been unable to find analytical expressions for 
the critical external field at which the macroion complex 
becomes unstable and had to resort to a numerical solu- 
tion of Eq. (|l|). The numerical procedure, however, is 
basically identical to the one without an applied exter- 
nal field, hence numerically very inexpensive. The result 
is depicted in Fig. ^. The roughness of the curves is 
not a sign of a numerical problem but rather due to the 
physics of the problem. Depending on the precise value 
for N, the geometrical arrangement is more or less sym- 
metrical, resulting in large changes in the dipole moment 
when TV is changed by only 1. Without an external field, 
this dipole moment is not relevant, and all quantities are 
smooth functions of N. This is no longer the case now. 

In Sec. IV we have shown that it is always possible to 
bind at least 2\Q\/q microions in a locally stable manner. 
Thus, for N < 2\Q\/q a finite electric field is necessary to 
break up the complex. For larger N, however, the critical 
field may vanish, explaining the division of the diagram 
into two separate regions by the line N — 2\Q\/q. 
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FIG. 5: Electric field E above which a state with N microions 
bound to the macroion is destroyed. (All labels in the 

figure are to be understood as \Q\/q.) 



VII. CHARGE DISTRIBUTION BETWEEN 
TWO MACROIONS 

Next, we apply our approach to the case of two well- 
separated macroions, with charges Q\ and Q2, and radii 
Ri and i? 2 , respectively, together with M = \Qi + Qi\/q 
microions so as to achieve charge neutrality. According 
to Fig. U and Eq. (|l^), there are many different pos- 
sibilities for distributing those M particles among the 
two macroions in a locally stable way. The naive result 
is the one where each dressed macroion becomes charge 
neutral. However, it is possible that 8 microions will 
be transferred from the first macroion to the second (if 
8 < 0, \6\ ions are transferred in the opposite direction). 
Such a state is called "ionised" @ Hi- Neglect- 
ing interactions among the two macroions, since they are 
well separated, and applying Eq. (||) gives the potential 
energy V(S) 



V(S) = [±f (\Qi\/q + S) + ^-J (IQ2I/3 - 5)} 
+ ^ e [(\Qi\/q + 5)j± + (\Q2\/q~S)^} . (17) 

We assume that the ions are at very low but finite tem- 
perature such that the system can break out of a locally 
stable state, and to find its lowest energy state ]35| . This 
state is easily computed from Eq. ( |l7|) and provides a re- 
lation between Qi, Q2 and S for the ground state. This 
result is most easily displayed when Q2 is expressed in 
terms of the other parameters, 



IQ2I =5+(S- 
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1 + S 
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R2 



(18) 



This curve is depicted in Fig. |6| for different ratios of R\ 
and i?2- For two identical macroions we find that the 
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FIG. 6: Overcharging (> 0) or undercharging (< 0) of the first macroion (in units of q) as a function of the charges Qi and Qi 
of the two macroions. The radii of the two macroions are R\ — R2 (left), Ri — I.5R2 (centre) and Ri = 2R2 (left). The lines 
are for the lowest-energy configuration. 



lowest-energy state is the non-ionised one whereas for 
Qi 7^ Qi but R± — i?2 the macroion with higher charge 
attracts more microions than naively expected. This is in 
agreement with earlier results for R\ = R2 where the cor- 
relation effects in a three-dimensional layer were approxi- 
mated by the effects in a two-dimensional Wigner crystal 
and confirmed by numerical simulations ]l6|, [n], |l9) . For 
R\ i= R2 the macroions are uncharged only if 

9i = 9i (19) 

Other configurations than the lowest energy state can be 
excited thermally, and, due to local stability, can persist 
for relatively long times. 

In a strict mathematical sense, at finite temperature 
microions cannot be bound to a three-dimensional struc- 
ture like a sphere. This is in contrast to a rod or a planar 
geometry |t| . If the temperature is not too high, or equiv- 
alently, if the electrostatic coupling is strong enough (i, e., 
large \Qi\, IQ2I and q as well as small R), the methods 
presented in this paper can still be applied in an approx- 
imate fashion. Microions stay very close to one macroion 
for most of the time, before they hop (i. e., move within 
a time that is short compared to the time that they re- 
main effectively bound) to the other macroion. If we 
neglect the short hopping phases, the probability for a 
given ionisation level S is then given by the Boltzmann 
factor, P(A) oc exp[— (3V{5)\, and is easily evaluated nu- 
merically for arbitrary parameters. 



VIII. CONCLUSIONS 

To conclude, using general stability concepts we have 
shown that for a charged macroion a stable overcharged 
state persists at zero temperature. The number of glob- 
ally stable microions is given by Eq. (|^), the number of 
locally stable microions is given by Eq. ([i2|). Physically, 
global stability means that a random arrangement of mi- 
croions will move to form a layer around the macroion, 
whereas local stability means that a layer that exists due 
to initial conditions will persist forever. An applied ex- 
ternal field, as used in electrophoresis, decreases the sta- 
bility of the microion cloud, as demonstrated in Fig. |^. 
We have also applied this approach to a system of two 
macroions. In its ground state each of the two macroions 
acquires a nonvanishing net charge, see Eq. ([l8]), unless 
Eq. @ is fulfilled. 

The value of our approach lies in the fact that it is ex- 
act, and by that clarifies the effects caused by finite size, 
finite charge, correlations and electrostatic interactions. 
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At very large number of microions, the curvature of the 
macroion surface seen by each microion becomes negli- 
gible. Then and only then the microion layer becomes 
two- dimensional. 

Imagine the bare macroion moving against the direction 
of the field to infinity, and the microions in the direction 
of the field. 

This sign is easily understood by noting that the ener- 
gies due to macroion-microion interaction as well as due 
to microion-microion interaction become smaller in mag- 
nitude when the fc-th microion is moved away but the 
energy due to the external field becomes larger in mag- 
nitude. 

Note that this state is not globally stable even though 
it is the state of lowest energy. In contrast to Sec. [II 
where no energy barriers existed, they do now 



